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^ ■ Abstract 

The one-loop fermionic contribution to the probability of an instanton transition with fermion 
number violation is calculated in the chiral Abelian Higgs model in 1+1 dimensions, where the 
fermions have a Yukawa coupling to the scalar field. The dependence of the determinant on 
fermionic, scalar and vector mass is determined. We show in detail how to renormalize the fermionic 
, determinant in partial wave analysis, which is convenient for computations. 
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I. INTRODUCTION 



The interest in the chiral Abehan Higgs model in 1+1 dimensions hes in the fact that it 
shares some properties with the electroweak theory, but is much simpler and may serve as 
a toy model. One of the most interesting common features of the two gauge theories is the 
fermionic number non-conservation m. Both give rise to instanton transitions, leading to 

nn 

the creation of a net fermion number due to an anomaly 12, Both theories contain finite 

n □ 

temperature sphaleron transitions |^ - |9||. 

At zero temperature, zero fermionic chemical potentials and for a small number of par- 
ticles participating in the reaction, the probability of the process can be computed using 
semi-classical methods. In general, the result is a product of the exponential of the classical 
action e'^"' and the fluctuation determinants. The latter factor includes the small pertur- 
bations of the fields around the instanton configuration, and in many cases may only be 
computed numerically. 

Quite a number of computations of determinants in 1+1 dimensions can be found in the 
literature^. In particular, the determinants have been calculated for the vector and scalar 
«e.d «uct.a.o,. around the .n.tanto. in Q. as well a. fo. the fe^.onie ones ,n M. where 
it was assumed that fermions have no mass term and no interaction with the Higgs field. 
However, to our best knowledge, no computations incorporating the Yukawa coupling of 
the fermions to scalar field have been done till now, neither for realistic case of electroweak 
theory nor for the chiral Abelian Higgs model ^. 

The aim of the present work is to partially fill this gap, calculating the fermionic deter- 
minant in the 1+1 dimensional case, where the fermions interact with the Higgs field in a 
similar way as in the electroweak theory ^. 

This calculation is somewhat delicate because of the difficulties occurring in regularization 
and renormalization of chiral gauge models beyond perturbation theory. Furthermore, an 
analytic solution to this problem cannot be obtained, since even the classical instanton 



^ For 3+1 dimensional computation without Yukawa couplings see the seminal paper by 't Hooft _ 
^ The determinants in the high temperature sphaleron transition in 1+1 dimensions were computed in |3,|9|. 
^ Similar studies have recently been performed for other models. In a sup ersymmetric theory in 2+1 



, supe 



dimensions the calculation is simplified by a supersymmetric constraint |13|. The fermionic contribution 
to the vortex mass has been calculated in a model resembling the (non chiral) Abelian Higgs gauge theory, 
where the fermion couples to the absolute value of the scalar field . 
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profile, given by the Nielsen-Olesen string solution Q|, is not known analytically, apart 

from the special case where the Higgs mass equals the vector field mass [l^. Nevertheless, 

we use analytical methods as long as possible before moving on to numerical computation. 

r~~i 

We will use a numerical method developed in ^m], extended to our case. 

The paper is organized as follows. In section 2, the model and its basic features, such as its 
vacuum structure, anomaly, instanton configuration and fermionic zero modes are discussed. 
In section 3, we study and compare the 1-loop divergences occurring in this model in various 
regularization schemes. In section 4, the method of |l6| to calculate determinants is discussed 
and applied to our case. In section 5 we present in some detail the numerical procedures and 
give the results of the determinant computation. Finally, conclusions are given in section 6. 

II. THE MODEL 

The model we consider here contains a complex scalar field with vacuum expectation 
value f ; a vector field A^, and Uf fermions , j = 1, rij: 

-V{<f)) + i + - ^/^r l^vl>^0. (1) 

The charges of the left and right-handed fermions differ by a sign, = — = ^ and the 



symmetry breaking potential is chosen to be V[(j)] = ^(l^l — v ) . In the following, we use 
the Majorana representation for the 7-matrices: 

n i \ 

, 75 = 7o7i- (2) 





^0 -^\ 
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Note that we do not reduce generality in considering Yukawa interaction between identi- 
cal fermions only^. In principle another mass term of the form M\E'^7o\I' + h. c. could be 
added to the Lagangian (Q). It is compatible with gauge and Lorentz invariance but breaks 
fermion number explicitly. As we are interested in instanton mediated fermion number 
non-conservation, we will not consider this term. 



A more general interaction could be written in the form + h.c. but the matrix can 

always be diagonalized and made real trough redefinition of the fields "if^ . 
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This model has been studied as a toy model for the fermionic number non-conservation 
in the e.ect..owea. .heo.y in a „un,be. of pape.. see. e., Q-Q. Q. 

The particle spectrum consists of a Higgs field with mass = 2Af ^, a vector boson of 
mass rriw = ef, and rij Dirac fermions acquiring a mass = f^v via Yukawa coupling. 
The model is free from gauge anomaly. There is, however, a chiral anomaly leading to the 
non-conservation of the fermionic current, 

Uf nf Uf 

j=i j=i j=i 

with a divergence given by 



The vacuum structure of this model is non-trivial Taking the Ao = gauge and 

putting the theory in a spatial box of length L with periodic boundary conditions, one 
finds that there is an infinity of degenerate vacuum states |n), n G Z with the gauge- Higgs 
configurations given by 

27rn 

Ai = —-, (P = ve'—. (4 
eL 

The transition between two neighboring vacua, described by an instanton, leads to the 
non-conservation of fermion number hy nj units. In this paper we consider rij to be even. 
The case of odd nj, resulting in the creation of an odd number of fermions, is analyzed in 



2l|. 



A. Lagrangian in Euclidean space 

As the tunneling is best described in Euclidean space-time, we review here the corre- 
sponding equations and conventions. 

The Lagrangian may be rewritten in Euclidean space: 

£^ = -J^,F,, + ^W^''^{^,-^^-^,A,)^>^ + V{ct>) 

+ ^(D^0)t(D,0) - ^pW^^^^* + z/^f1^v^^0, (5) 

with = d^ — ieA^, •Jq = i'jQ and 7f = 71. The fields \1/ and \1/ are independent variables, 
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(6) 



and the gauge transformation reads: 
For comparison, the Lorentz transformation is: 

^(x) ^ ^'(x') =^a;^(a-V), 

with As = exp(i7^|) being the rotation matrix in two dimensions. 
B. Instanton 

The instanton which describes the tunnehng between the states |0) and |?7,) is simply 
the Nielsen- Olesen vortex with winding number n which is a solution of the Euclidean 
equations of motion in two dimensions. In polar coordinates (r, 9), the field configuration 
reads: 

(j){r,e) = e^"VW = e^"V(^), (7) 
A\r,e) = e'^PA{r), (8) 

where r is the unit vector r = (cos^, sin^) and e^^ the completely antisymmetric tensor 
with = 1. The functions A and / have to satisfy the following limits: 

/(r) ^ crl'l, 
fir) ™ 1, 

Air) ^ 0, (9) 

A(r) --. 

er 

Passing to dimensionless variables 

TTl ~ Tfl ~ T 

A = — A, (p = — (j), r = — with m = v Xv'^ (10) 
e em 

reduces the number of free parameters. The equations for A, cj) are : 

- df {^dffAir)^ + 02 (^A{r) - =0, 

-^df (f9,0(r)) + A{r-)\ - 1 + ^i^irf ) 0(r) = 0, (11) 



with jj, 



The classical action is given by: 



2m; 




(12) 



The number AA^ of fermions created in the instanton transition can be computed by inte- 
grating (jnj over the Euclidean space: 



where q = J cPx^e^^F^^, is the winding number of the gauge field configuration. For the 
instanton configuration ()7p8|l . we have q = n. 

C. Fermionic zero modes 



According to the index theorem (see for example '26|), the Dirac operator in the back- 
ground of the instanton satisfies the following relation: dimker[i^] — dimker[i^"''] = n. As 
the instanton in 1+1 dimensions coincides with the vortex, these zero modes may be found 

n 

by carrying out a similar analysis as in j23| ; where the fermionic zero modes on the Nielsen- 
Olesen string were analyzed for non chiral fermions. In this subsection we present the 
corresponding equations. 

The Lagrangian for the fermion j in the background of the scalar and vector fields may 
be written as £f, = K^"^^ , where 



In the following, the family dependent Yukawa coupling will be replaced by / keeping in 
mind that there is no mixing between different fermionic generations. 

The zero modes are the regular normalizable solutions of the equation K'^ = 0, with 
Afj_ and given by ()7l8p ^. Using polar coordinates and performing the substitution \E' = 

^ However, in the massless case (/^ = 0), a logarithmically diver gent wave function is generally kept as a 
relevant solution. The reason is that its classical action is finite 1201 . 




(13) 





(14) 
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exp ^i^dp \1/ we get: 

-iFf(r)e-'''^ ipio (d_ ^ id^^\ ^ 

^ = 0, (15) 



where F = fv is the fermion mass. With the use of the phase decomposition \1/ = 
Sto=-oo e*™'^\l/™, equation (fT3j) can be rewritten as 

Ff(r)*T - (I - ^^^^^) = 0, 

— + -U^- Ff{r)^^—' = 0. (16) 



In our case, the analysis of j22] shows that for a vortex with topological number n < 
there are exactly |n| fermionic zero modes in the spectrum of K with m in the interval 
m G {—n + 1, .., 1, 0} and none in the spectrum of K"^ . For n > there are no zero modes 
in the spectrum of K, but n in the spectrum of K"^ . 

For the case of n = — 1 studied below the explicit form of the zero mode is given by 

^0 (r) = ^o^(r) oc exp |F/(r') + |A(r')} rfr'^ • (17) 



Note that for massless fermions {F = 0), the zero mode decreases as for large r. It i 



IS 



therefore not normalizable and has a divergent action. This behavior differ form the case 
of H. Q. a. where „,a..es. fe^ion. of charge, e were co,.,dered. I„ the. ease the 
fermionic zero mode decreases as - for large r and has a finite action. 



D. Determinant 



Due to the presence of the fermionic zero modes, instanton transitions imply the creation 
of a net number of fermions. In the following, we will be interested in the creation of one of 
each type of fermion, for which an instanton of charge n = — 1 is needed. The corresponding 
transition probability is proportional to det'i^", where the prime means omission of the zero 
eigenvalue in the calculation of the determinant. 

It is well known Q that the eigenvalue problem for the operator K is ill defined. Con- 
sequently one has to consider the Laplacian type operators K'^K or KK'' which have the 
same set of eigenvalues (except for the zero modes). Then det'i^ is defined up to a phase as 

7 



det'[i^'] = det'lK^KY^'^. The explicit expression for the operator K^K reads: 

K^K = (18) 

if[ct>*{eA, - leA^) - {ido + P\ct>{r)\' - {8^ + A^)^ + f e^,9^A, 

The fermionic equations of motion, for instance equation (fT3j) . remain unchanged after the 
variable changes (fTn|l . if / is replaced by //e and e is set to 1. The only free parameter in 
the bosonic sector (fT^ is /i, while there is a second parameter in the fermionic sector: the 
Yukawa coupling /. 

In conclusion, we are left with two dimensionless parameters, and the determinant can 
be calculated function of 

^ = ^^V2 and JL = = f (19) 

mw ruH e uih v2A 

Obviously the determinant, being a product of an infinite number of eigenvalues, is a 
divergent quantity. In the next section we discuss its regularization and renormalization. 

III. REGULARIZATION AND RENORMALIZATION 

For perturbative calculations, the dimensional regularization is best suited. However, as 
has been observed in it is not applicable to the computation of the fermionic determi- 
nant because the continuation of the instanton fields to a space with fractional number of 
dimensions is not uniquely defined. Nevertheless, we discuss the dimensional regularization 
to fix the meaning of the Lagrangian parameters in section 3.1. In section 3.2, we consider 
another regularization scheme based on partial waves decomposition. It permits to exploit 
the spherical symmetry and turns out to be convenient for numerical purposes. In appendix 
A, we consider the Pauli-Villars regularization used in jl| and prove its equivalence with the 
partial waves procedure. 

A. Dimensional regularization 

The Lagrangian depends on four parameters; the charge e, the scalar coupling A, the 
scalar mass rriH, and the Yukawa coupling /. The model under consideration is super- 
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renormalizable. In order to use the dimensional regularization we have to define the 7- 
matrices for an arbitrary number d — 2 — e oi dimensions: 

tr(7V) = 2^'^^ iy = 0,l,...,d-l. (20) 
The definition of the 7^ matrix is ambiguous, we follow here the usual definition: 

[7^71 = 0, u^2,...,d-l. (21) 

The physical parameters e, A, mn, f in two dimensions are related to the d-dimensional 
parameters e^, A^, rnnd, fd by: 

ed^eii^~^, Xd^Xfj,'^''^, mHd^rnHlJ^~K fd^flJ^~^- 

We will work in the gauge. The complex field is written a& (f) — v -\- h -\- iip, where h 
and </? are real. The gauge fixing term is 



where 



G[A,h,<p]^^{d^A^-^ev<p). 



The Lagrangian for ghost fields c is 

^ghost — ^ 



1 



(22) 



(23) 



(24) 



In the following we will work in the minimal subtraction scheme. The only divergent pa- 
rameter is the Higgs mass uih- A straightforward computation gives the relevant part of 
the effective action, 

S^nt - I d'xl (0^ - v') 6m' (25) 

with 



Sm' — — 
An 



MS 



(26) 



where = ^ ^ 7 + ln(47r). In the minimal subtraction scheme we subtract the coun- 

terterm 

with Sm^^ containing all terms in (j2Sl) proportional to 

For the photon propagator, the bosonic loops do not introduce any renormalization. 
However, as is well known Q], there is a finite contribution coming from fermionic loops. 
Because of the ambiguities in the definition of 75, dimensional regularization breaks the 
chiral gauge invariance and a term ^^4^ needs to be added to the action. The complete 
counterterm action to be subtracted from the initial action (jS)) reads: 

Scount ={ld'x (^-^Al + ^ (0^ - v') 5m^^ } . (27) 
B. Partial wave regularization 

The spherical symmetry of the instanton suggests that partial wave expansion can be 
used. The eigenvalue problem decouples into one-dimensional differential equations. In this 
section we discuss a natural way to regularize the partial waves. We consider here only the 
fermionic sector. 



1. Partial wave expansion 



We may write det[K^K] path integral: 



det[K^K] = j Vr]Vf]exp 



(f'x rfK'^ Kr] 



Partial wave decomposition is defined as follows: 

00 00 



(28) 



(29) 



The regularization is done by putting our system in a finite spherical box of radius i?, and 
cutting the sum over the partial waves at some m = L. After performing the partial wave 
decomposition, the regularized action reads: 



S{R,L)= / 2nr dr ^ f]"'{r)M"'^{r)ri^ 



(30) 
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with 



27r 



(31) 



From the general expression (|T9|l for K'^K, we get for the vacuum: 



K^K 



- ^2 - df 

- ^2 - df 



(32) 



After phase decomposition we obtain a diagonal matrix in both spinor and partial wave 
space: 



ml jrml 



5™' 11, 



1 9 m 



where 112 is the identity in spinor space. The radial eigenvalue equation in vacuum reads: 



(33) 



n/rrnm m _ \2 m 
^^^vac '/A — ^ '/A ' 



(34) 



with boundary conditions 



ry™(0) = 7]"'{R) = 0, m ^ 0, 
r/°(0) = l, r/°(i?) = 0. 

From the relations ()33j) - (j!?H|) the free propagator may be derived 



(35) 



1 



A2 



/m(Fr) 

Im(FR) 



[K^{Fr')I^{FR) - Im{Fr')K^{FR)] , r < r' , 



(36) 



2vr 1 MJrO [K^^Fr)Im{FR) - lUFr)Km{FR)] , 
It allows us to treat the interaction terms present in ()19j) by standard diagrammatic methods. 



2. One-loop divergences in partial waves 

As we have already seen, the fermionic parameter / needs no renormalization. However, 
the mass of the scalar Higgs receives divergent contributions from fermionic diagrams. The 
partial wave regularization can't be introduced at the level of the fermionic Lagrangian 
(0), but only at the level of the squared determinant (j2Hl)- One does not expect that the 
counterterms derived from the initial Lagrangian are sufficient to remove all infinities in 
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). Hence, we recalculate the counterterm action (see appendix B for details) needed to 
renormalize (PHj) . The result is: 

L 

m=--L 

= J2 ^27rrtT[G-{r,r)]ff{\<j>\'-v') + ^e,^^,A^dr, (37) 

where the S^^^^ are finite for each m and only the sum is divergent in the limit L —>■ oo. Note 
that the counterterm (jH7|) is non-local. This is due to the non-locality of the partial wave 
regularization procedure and may be checked to be correct by comparison to Pauli-Villars 
regularization, see appendix A. 6. 

For small constant background fields, ^7\i leads to 



S^Lt = [ {f (l0(r)P - v^) + le,.d,A^{r)) d^r^ [l + log 



(38) 



In order to get results in the MS'-scheme from those calculated in the partial waves, we 
calculate the difference 5S^^^^ between the effective action found in these two schemes^. 
The result reads 

, ^ (log (i^) - (1) J / - + (39) 

where 

S^f = log (^) / ^^^e.AM^)- (40) 
In comparison with the dimensional regularization, a supplementary divergent term in- 
volving gauge fields Sgf has appeared. It also arises when using Pauli-Villars scheme (see 
appendix A. 3) and is an extra divergence of the action ()28p in comparison to the initial 
action (0). If we Wick-rotate Sgf back to our initial Lagrangian in Minkowski space-time, 
it gets an extra factor of i, the action becomes non-hermitian and breaks unitarity. Because 
of this, Sgf must be subtracted completely. 

For the photon propagator, as in dimensional regularization, we have to subtract from 
the effective action the term 



e2 rR 



SZntiR) = ^l d\A\r) (41) 



to recover chiral gauge invariance (see appendix B.l). 



^ The counterterms found in the dimensional regularization have to be multiplied by a factor of 2 because 
we are dealing here with the squared operator K'^ K. 
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3. Regularization and renormalization in partial waves 



From the counterterm ()38|). we see that the initial theory is recovered in the hmit 
oo. The summation over partial waves and the limit L —>■ oo has to be performed first and 
the infinite volume limit must be taken only after having removed the infrared counterterm. 

The explicit expression for the counterterms in the case of space dependent background 
is obtained in integrating (jHTj) and the renormalized fermionic determinant may formally be 
written as: 

detren[K^K] = (42) 



lim lim 

i?— >oo \ L—>oo 



exp{-SZntiR)} 



This prescription differs from the one of where the limit i? — oo is taken first. It is 
shown in Appendix C that the order of the limits is crucial. 



IV. DETERMINANT CALCULATION 

After the partial wave decomposition (EHHSU), K^K{r,e) was expressed in terms of 
M'™'(r). For our purposes, the case where M''" is diagonal in partial wave space (M'*" = 

S^M"^) is sufficient^. The determinant may be calculated as: 

+00 

det[K^K] = Yl detlM""]. (43) 

m=— oo 

We are left with the much easier problem of finding the determinant of one-dimensional 
operators, which may be addressed with the following theorem Q]: Let us consider two 
operators Oj = —d^ + Wi{x), i = 1,2 defined in an interval of length R. Let i = 1, 2 be 
the solution of Oi'^i = with the boundary conditions 

*,(0) = 0,^:(0) = 1, z = l,2, (44) 

we have: 



det 



Oo 



^ We are mainly interested in the case n = — 1, where one of each type of fermions is created. In this 
particular case, the operator Mj™g( is diagonal in partial wave space. Note that this point is not crucial, 
as explained in section 4.1, the determinant may be calculated in the non-diagonal case as well. 
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A. Treatment of radial operators 

We follow here the method developed in to calculate determinants. Note that here 
we will first consider the radial problem for to i? where 1 and the limit R oo will 
be taken afterward. 

In the present case, even if M™' is diagonal in partial wave space, it is not diagonal in 
spinor space. The theorem ()45|) needs generalization to two coupled second order differential 
equations. We are interested in the ratio between the operator 



in the instanton background and the vacuum operator M'^'^"''^, which is assumed to be 
diagonal. Let us define the matrix ip'^ = 1,2) and 'ip^'^^ as the solutions of the 
following differential systems: 



(46) 



with boundary conditions 



lim 




= 1, 


lim 




1 vac,in 


1 vac,m, 


lim 




= 0, 


lim 




/ vac,m, 


1 vac,m 
^R 



The determinant is then given by: 

det[M^] det[^,7(i?)] 



0, 

1. (47) 

(48) 



det[M™'--] 

The remaining determinant is just the usual determinant for 2x2 matrices. It is an easy 
exercise to reproduce step by step the demonstration of 2^ in this more general case. 

The vacuum operator which is given in ()33|) . has an analytic solution ^t'^''"'^'^ = I^(^Fr). 
For the instanton {n = —1) configuration we get: 



= M^,=F (-fir) - -f{r) + eA{r)f{r) 
\ r 

1 <9 w? ^9,9, X e ^ , ,9, , Air) ,^ , 

= - + ^ + F^f\r) + -e,^d,A^ + -A\r) - me^. (49) 
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The solution ^^(r) needs to be computed numerically. To this aim, it is convenient to make 
the following substitution: 



(50) 



The determinant is then evaluated with 

det[M™] 



(51) 



det[M^"^'™ 

In terms of the functions /i^, the equation takes the form of an ordinary quantum 
mechanical equation with potential Vij{r): 



I n , MFr) \ d_ 
dr"^ \r Im{Fr) J dr 



hT,{r) = V,,{r) (4^. + /,™.(r)). 



(52) 



The effective potential Vij{r) in the background of the instanton is given by the following 
expressions: 



2 . .2. s N ^{r) 2^^(^) ^{r) Air) 
l^n(r) = {f{r) - l) - + e^—^ - e-^ - me^. 



Vi2{r) = V2i{r)=eFf{r)A{r)-Ff{r) 



Ff{r) 



(53) 



V22{r) = F^ {f{r) - 1 - + e^— ^ - + me^^. 

2r 4 2 r 



The functions h^{r) can easily be found numerically from with the boundary conditions 



(54) 



h,,iO) = 0, t,j = l,2, 
/i^(0) = 0, t,j = l,2. 

For m = we have to remove the zero- mode present in M^^^^. In this case, it is possible 
to diagonalize the operator M™^^ with the substitution = \E'i ± \E'/j: 



1 d 



+ F'f\r) + -e,,d^A, + -A\r 



±i^(-/'(r)-^/(r)+eA(r)/(r)^ 



^± = M?^± = 0. 



(55) 



The fermionic zero-mode is contained in M^. We calculate det[M^] with and det'[M^] 
as in [iG^ : 



/ r 71 rG,inst-\ 



det'[M 



^det[Mj^"^* + A^]U.=o d 



detiM"'™'^] 



det[M 



h\R). 



(56) 
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In the last relation h^{r) is defined through = vl/^'^'=(l + with \l/^_^*(r) being a 

solution of 

with the boundary condition (j44p and '^'""■'^[r) = I^^Fr) being the solution of 



B. Ultraviolet divergences 

A possible way to calculate the counterterms S^^^^{R) is given in (jHTj). We need to 
integrate numerically the Green's function multiplied by the potential U{r) = — 
f ^) + ^Sfj^iyd^A^. As the Green's function is not smooth, for numerical calculations, it is more 
convenient to solve the related differential equation^ 



+ - + 2 



dr"^ \r Lfn{Fr) J dr 



STountir) = 2U{r), (57) 



with the boundary conditions 5*^.(0) = 0, S'j^f{0) = 0. For the instanton configuration we 



have: 



U{r) F\Pir) - 1) - f (A'(r) + ^). (58) 



V. NUMERICAL PROCEDURES 



In this section we describe the numerical methods used in this work. First the background, 
namely the well known Nielsen-Olesen vortex is considered. The method used here to find 
the profile is explained briefiy. In the second part, the calculations related to the fermionic 
determinant are discussed, namely the integration of the differential equations, asymptotic 
solutions, subtraction of divergences and treatment of zero-modes. The renormalization and 
convergence of the different limits are checked and finally, results for the determinant are 
given. 

A complete derivation is found in the extra factor of 2 in front of U{r) conies from the trace in spinor 
space. 
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FIG. 1: Instanton profile for rriH/mv/ = 2, classical fields (j){r) (dashed), A{r) and their asymptotic 
forms (thin lines). 

A. Background 

The instanton profile may be found with a shooting method (see for instance 
boundary conditions at r = are of the form: 

v4(0) = 0, v4'(0) = 6, (/)(0) = 0, (^'(0) = (59) 

where the parameters h and (3 are found imposing the limits (jHl) and ()10|). We start the 
numerical integration at r ~ 10"'' instead of r = 0, where some trivial divergences occur, 
and use a small r expansion for and A: 

A{r) = br- —r^ + 0{r^), 
8 

0(r) = pr- ^ii±^r3 + 0{r'), (60) 
8 

valid for n = —1. The numerical integration is done with 32 decimals, and to get an 
accurate^ profile, the boundary conditions have to be specified within an accuracy of order 

^ The accuracy can be checlced by calculating the instanton number or the action of the instanton for 
= 1 that is known to be ttv'^ 0|. The results of the numerical integration agrees to 13 decimals with 
the action in the latter case and at least 7 for the instanton number in any case (see figure ^ and table 

rvTi . 
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B. Fermionic determinant 



For the fermionic determinant, the task is to solve equations ()52|) and (j37|) . These equa- 
tions are completely symmetrical under the change m —>■ —m therefore only positive m need 
to be considered. As the solution to the fermionic equations in the asymptotic instanton 
fields is known, it is sufficient to integrate numerically to r ~ 15 and glue the asymptotic 
solution 

^^(r) = A/^_i/2(Fr) + 5ir„„i/2(Fr), 

^^(r) = C/„+i/2(Fr) + DKm+i/2{Fr). (61) 

The constants A, B, C, D are determined in imposing the continuity of \l/™(r) and its first 
derivative. The numerical integration, like for the vortex, starts at e ~ 10^®, where the 
boundary conditions are found by calculating the power expansion for the K^: 



Having found the /i™, we calculate the partial determinants with (pSj) and subtract to each 
wave the partial counterterm S'^^^^ found with (fHTj) as prescribed in (jl^ . 

Numerically we store the value of the determinant for ~ 50 different values of system 
radius i = 1, ...,50. After renormalization, the partial determinants det[M'"] have to 
decrease at least as so that the product over m remains finite. This is checked in figure 
2. Using this property, for each Ri, we calculate the partial determinants from m = 1 to 
m = L ~ 30 and fit them with an inverse power law: 

A f r Mm] _ const2 const3 const4 

m-^ 

This approximate expression is then used for m = L ~ 30 to infinity. 

This completes the limit L ^ oo and we may consider the limit R ^ oo. At this point 
the determinant still depends on R (see figure 3), according to ()42|1 . the infrared counterterm 
1)4111 have to be subtracted. The renormalized determinant becomes approximately constant 
for typically 10 < FR < 100 and FR can be chosen in this range. Keeping in mind that for 
large FR higher partial wave should be considered, it is expected that the result becomes 
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FIG. 2: Logarithm of the partial determinant log(det[-/Vf™]) as a function of angular momentum 

2 

m, for the parameter values = 1, = 0.1. The upper small dots shows the results for 
determinant before renormalization, and the lower ones after. The lines show respectively 1/m 
and l/m? behavior. 
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FIG. 3: Logarithm of the determinant log(det[M]) as function of the system radius R, before 



(dashed) and after subtracting the infrared counterterm (|1T]) 



1, 



rriH 



0.1) 



inaccurate at large FR. Fortunately the determinant converges very fast as FR oo and 
is found to be constant up to 4 decimals for typically 20 < FR < 40 from were the result is 
extracted. 

For m = the zero-mode in M? has to be removed in the determinant calculation. This 
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is done with()56p. where the derivative is approximated as 

= 4r2h\R) ^ (63) 
oA A 

To get an accurate result, we take A^ of the order (10~^)F^ and perform the computation 
of h^{R) for some (~ 10) different values of A. These results are fitted to extrapolate the 
value of at A = 0. 



C. Results 

We first note that det'[M] has dimension of mass^^ from ^ in (jKHjl . The fermion mass F 
may be used to obtain a dimensionless quantity F^^det[M]. The results for F ^/det/[K^lC\ 
are plotted in figure |3] The logarithm of the partial determinant log(det[M™]) behaves as 
— and after renormalization as see figure El It becomes constant at large R after 
subtraction of the infrared counterterm, see figure |H1 

The behavior of the determinant for small fermion mass F is a power law, see figure El 
This comes from the partial determinant det'(M^) where we remove the zero mode and can 
be checked with some analytical approximation (see appendix D). The accuracy of the value 
for the determinant is estimated to be of the order 10^'^ but may be less for < 10~^. 



20 



- - - 1/5 
0.1 



-2.5 -2 -1.5 -1 -0.5 

Log [F/mn] 




-2.5 -2 -1.5 -1 -0.5 

Log [F/mn] 



FIG. 4: Logarithm of the determinant log [F \J det' [iiTt ] ] as a function of the dimensionless fermion 

2 

mass (horizontal axis) for different values of = -^r-- The different values of the determinant 
are fitted to few percents accuracy with the following expression: F ^/ det' [K^ K] = 1.62/x^/^°F^/^ + 
(-4.60 + 3.71/xV5 - 0.632 In /i) F+ (6.97 - 6.76/xV5 + 0.866 In /j) ln(l+F) in the interval < F < 3 
and 0.05 < < 8. 

VI. CONCLUSION AND OUTLOOK 

In this paper, wc have studied an instanton transition in the chiral Abelian Higgs model 
with fermion number violation and computed the fermionic determinant taking into account 
the Yukawa couplings. 

The dimensional regularization has been used to fix the meaning of the Lagrangian pa- 
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FIG. 5: Determinant Fy/det'[K^K] as a function of the dimensionless fermion mass for the 
case ^ = 1. The values are fitted with the power law: F^detMl = IMA (7^^)^^^- 





Scl 


b 


P 


1/10 


0.6388286986270 


2.557798983491183 


5.756251019029544 


1/5 


0.7259086109970 


1.554461598144364 


3.541849174468259 


1/3 


0.8008642959782 


1.081478368993385 


2.496453159112955 


1/2 


0.8679102902678 


0.812560321222651 


1.901012558603257 


2/3 


0.9199259150759 


0.663981767654766 


1.571374124589507 


1 


1.0000000000000 


0.499999999999919 


1.206575709162995 


2 


1.1567609413307 


0.308286653343485 


0.777359529040461 


4 


1.3405945494178 


0.189926436282935 


0.508674018585679 


6 


1.4612151896139 


0.142825844043109 


0.399789567459296 


8 


1.5526758357349 


0.116536242666195 


0.338046791533589 



TABLE I: Results for different — ^: classical action {Sd), logarithm of the fermionic determinant 
for dimensionless fermion coupling g = 0.1 (detir=o.i) and the boundary conditions at r = for 

the instanton profile (b and P). 
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rameters. The numerical calculations have been performed in the partial wave scheme, and 
the Pauli-Villars regularization is studied for completeness in Appendix A. 

In the limit of massless fermion (F — > 0), our results can't be compared to the calculation 

n n n 

of [18||, j20||. Fermions of electric charge equal to the scalar field charge e where consid- 
ered in these previous references, whereas we considered pairs of fermions with half-integer 
charge |. The instanton transition probability vanishes as F^^^ in our case whereas it is 
finite in the case of integer fermionic charges. As noted in section 2.3, there is no fermionic 
zero mode in our case if the fermion mass is set to zero. It is therefore not possible to create 
massless fermions with an instanton of charge n = —1. The fact that the probability to 
create fermions vanishes in the massless limit confirms this observation. 

As can be seen in (I13|) . considering only one family of fermions leads to the creation of one 
single fermion. This process seems to be possible in two dimensions although is it forbidden 



24| . This is an important feature of this 



in four dimensions because of the Witten anomaly 
model, which is addressed in j^lj. 
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APPENDIX A: APPENDIX: PAULI-VILLARS REGULARIZATION 

We compare here the Pauli-Villars regularization of Q to the MS regularization and 
partial wave regularization. The partial wave regularization shows some unusual features 
such as non-locality, see equation (jHTj) . and a renormalization of the gauge field action, see 
equations ()39M40|) . In order to understand better their origin, let us compare the partial wave 
and the well known Pauli-Villars procedure. In Pauli-Villars regularization, a determinant 
can be calculated as in ll|: 



dpt IK^K 1 det[irtK^,,] det[KtKo + M^] 

det.e.[K K^,,] = det[irtKo] det[Ktir^,, + MT ^ ^ 

where K^Ka^ci, and K'^Kq are respectively the fermionic operators (fT^ in the background 
of the fields {A, 0) and in the vacuum. In order to determine all necessary counterterms 
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in this regularization scheme, one can consider small perturbations around the vacuum. In 
principle, the instanton determinant under consideration may have been calculated within 
Pauli-Villars regularization. However, the partial wave analysis is technically simpler for 
numerical computations. 

1. Effective action in Pauli-Villars regularization 

The potentially divergent terms may be extracted in calculating the first and second 
order terms in the Taylor development of the logarithm of ()A1|) with respect to the fields. 
In the sections A. 2, A. 3, A. 4, we calculate all the relevant functional derivatives and find 
their contribution to the determinant. The result is the following effective action: 

SZnt = log(l^) / ^^{le,.d,AAx) + f{mx)\'-v^)}. (A2) 

Note that the Pauli-Villars regularization is gauge invariant , however as in the partial waves 
pnjl . a new divergent term proportional to e^adfj,A(j arises. 

To make the link between Pauli-Villars and dimensional regularization in the minimal 
subtraction scheme, we may calculate the difference SS^/J^^^ between the effective actions 
()A2|) and ()25|) . This provides us with a way to interpret the Pauli-Villars parameter M in 
terms of the parameter /i coming from dimensional regularization: 

«-.(.,M) ^ (.og(f ) - (i) J / '^fm.)?-.^-) 

+ log(^j J ^^e^ad^A^ix). (A3) 
The renormalized fermionic determinant detren[K'' K] may be written as: 

detren[K^K] = lim {detreAK^K] exp {-SZnt}) ■ (A4) 

Al — >oo 

The counterterms calculated above can be checked to be sufficient, in calculating determi- 
nants of configurations of A and that contains small perturbations around the vacuum. 
This can be done analytically, see appendix A. 5. 

2. Functional derivatives with respect to the scalar field 

We consider first the contributions that lead to the renormalization of the scalar field 
mass. The corresponding divergent terms can be found in calculating the first and second 
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derivatives of ()A1|) with respect to the scalar field. The first derivatives read: 



6(t){k) 



log [d.eireg[K'^ Ka,4\) 



A=0,(l}=v 



log {detreg[K''KA,^]) 



A=0, <j)=v 



(A5) 



and their contribution to the logarithm of the determinant is 

"^'^ Mk) - V6\k)) ^ log {detreg[K^KA,^]) 



+ h. c. 



A=0,(t)=v 



= / d^x {v{(j){x) -v)+ v{(f)\x) ~ ^) } ^ log ( ^ 

with 5"^ the two-dimensional Dirac delta function. The second derivative reads 

6 6 



(A6) 



log {detreg[K^K{</),A)]) 



A=0, ib=v 



= |-5(A: - g) log (^^j+0(M-2), 
which gives the following contribution to the logarithm of the determinant: 



^ / Sx\(t){x) 



f P lo£ 



M2 



The contribution of the first and second order derivatives can be added to give the second 
term in ()A2jl . which represent a renormalization of the Higgs mass. 

3. Functional derivatives with respect to the vector field 

The first derivative with respect to the vector field 

5 



5A,{k) 



log [deireg[K^ Ka,, 



A=0, ip=v 



does not vanish and gives a contribution to the determinant of the form: 



(fx 



^fiad^Acr{x) . 



(A7) 



Note that / ^^ficrd^Aa-^x) is just the topological charge. For small perturbations around 
the vacuum with usual boundary conditions (like infinite space and finite energy, see 
this integral is equal to zero. However this is not true in general and in the present case this 
integral is equal to —1. 
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4. Photon mass term with Pauli-Villars regularization 



The regularization procedure we used is gauge invariant. As this is not completely trivial, 
we shall now check it evaluating the one-loop corrections to the photon propagator. This 
can be calculated with the second derivative of ()A1|) with respect to A or by evaluating the 
corresponding Feynman diagrams. Surprisingly, the result, in the limit where the photon 
momentum q goes to zero reads: 

6 6 



lim 

g2-o 6A^{k) 6Ap{q) 



log {detreg[K^KA,4>]) 



in 



(A8) 



A=0, <t>=v 

That is we get a mass term for the photon. It is important to note that in chiral gauge 
theories regularized with a non-gauge invariant procedure, this is a common feature. But 
here, unlike for instance in dimensional regularization, the regulator term is gauge 
invariant, and such a problem should not arise. Indeed performing further calculations, we 
can see that every term of the scalar field covariant derivative |((9^ — 26^4^)0 p receives a finite 
contribution from the fermion loop, so that this vector field mass term can be absorbed in a 
gauge invariant expression. This confirms that Pauli-Villars regularization preserves chiral 
gauge invariance. 

In the remaining of the section, the calculation of the fermionic contribution to photon 
propagator is presented in more detail. Three diagrams are divergent or constant when the 
photon momentum goes to zero. We do not present the full calculation by second derivatives 
of the action but only these three main contributions. The first diagram is the 1-vertex loop 



with ^Af, interaction: 

4 i-^ 




e 

4 27r 



tr(]l)5 



1 



1 



p2 ^p2 p2^p2^ 



log K (9(Af-^ 



(A9) 



The second diagram is the 2-vertexes loop with —ie'-j^A^d^j^ interaction: 



''f(^y^ - 2! 

e2 5, 



'' £p f tr((75)2) zp^ i{p + q)i 



flU 

A 2^ 



{2n) 



p2 + F2 + M2 {p + qf + F'^ + M^ 
0{q' 



M=0 



M 



(AlO) 



The integration over p is done by standard techniques. These first two diagrams cancel 
each other to 0{M~~^) + 0{q^)] but the third one gives some constant contribution. Let us 
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consider the 2-vertexes loop with —iefcp'j^A^ interaction; is considered to be in vacuum 
configuration (p = v and fv = F: 




2! 



M=0 



M 



(All) 



This gives equation ()A8|) . which violates at first sight the chiral gauge invariance. Let us 
consider now other terms involving scalar and vector fields get such contributions, namely 
—ieA^{(f)*d^(f) — (pd^cp*). Let us consider the diagram with two vertexes /7^9^(3?(0) + 
^75^^(0)): 



7<!>-- 



2! 



|^ + 0(M-) + 0(g^), 



M=0 



M 



which gives a contribution 



1 



(A12) 



(A13) 



4:71 f 2 

to the effective action. The next diagram is the mixed one and contains one vertex 
/7^9^(3?(0) +^7553(0)) and one —iefcp'j^A^; the product of these vertexes gives two terms: 

We drop the second one, which is not part of the scalar covariant derivative, and which is 
gauge invariant (up to total derivative): 



2ief 



2! J (27r)2 + ^2 ^ JVf2 (p + g)2 ^ ^ 
ie_%n^qu+^5^ [1 + 0(M-2) + 0(g2 



M=0 



M 



(A14) 



27r 2f2 

with q the momentum of incoming scalar field and q' the out-coming one. This gives a 
contribution to the scalar-gauge effective action of the form: 

ie 



27rf • 



(A15) 



It is now possible to resume the terms dMl in a manifestly gauge invariant term 

j^\{d^ — ieA^)(f)\'^ to be added to the initial scalar covariant derivative ||(9^ — ieA^)(j)\'^ 
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and the photon acquires a mass 



m^'-^'=^/e2t;2 + |l. (A16) 



I — loop j Z Z 

This mass can be expressed with the dimensionless parameters ()T9j) as ^^^^^^jj;^^ — ~ V 2^ + i^' 
which does not depend on the fermion mass. Note that in the case of massless fermions, a 
similar phenomenon appears (Schwinger mechanism jl7|). 



5. Determinants of small fluctuations 

We are checking here if the counterterms mentioned before are sufficient to get a finite 
determinant. In order to be able to do it analytically we will only consider some small 
constant perturbation and calculate the ratio of the determinants in ()A1|) : First let us take 
(p = V + 6(j), A = 0, and note that 6F = f6(f): 

^ det [K^K{<f), A)] det [K^K.ac + M^] 

= log 



det [Km^ac] det [K^K{(f), A) + M^] _ 
det[l{-d^ - df + {F + 5FY)] det[l{-di - Of + + M^)] 



det[l{-dl - df + F2)] det[l{-di - df + (F + 6Fy + M^)] 

In momentum space, we can rewrite the last expression as 

d^k { ie + (F + 5Fff + F2 + M2)2 

log 



(27r2) ^ V + F^f (k^ + {F + 6Ff + M^f 

which can be easily calculated to give 

As the logarithm of the determinant is the sum of all one loop diagrams, we have to make 
subtractions at this level. Clearly the second term of the counterterm ()A2|) removes the 
divergence of this determinant. Then if we take cf) = A^ = 5A^ small constant pertur- 
bations; it is easy to perform the same calculations to see that no divergent term occurs. 
Similarly if we take simultaneously = f + 50 and A^ = the calculation is more compli- 
cated but we recover once again the previous divergences. However we can see that, taking 
a specific configuration where e^^dfj^A^, is constant, (p = v and A = 0, we find a divergent 
contribution of the form: 

^^e^^a^A^ log I ^2 



which is subtracted exactly by the counterterm ()A7|) 
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6. Equivalence between Pauli-Villars and partial wave 
counterterms 



For small constant background fields we may compare the partial wave counterterm (|37p 
and the Pauli-Villars one ()A2|) The difference SS^^^^ between them relates the different 



cutoffs M and ^: 



For any background that approach vacuum at infinity, it can be shown that the Pauli-Villars 
counterterms are equivalent to the partial wave ones. We introduce a Pauli-Villars regulator 
in the partial wave counterterm ()H7j) : 



V f 2'Krti[G'^{r,r)]h{r)dr 

L i-R 

lim y / 2nr (tr[G^(r, r)] - tr[G'^(r, r)]) /i(r)cir 
J~r Jo 



in=—L ' 



with h{r) = + ^s^ydfj^Ay) and Gp the Green's function for a particle of mass F given 

in equation The sum over m is now convergent and we can take L ^ oo. The sum of 
the Green's functions reads: 

oo ,00 

V (tr[G'^(r,r)]-tr[G'^,(r,r)]) = - V {Im{Fr)K^{Fr) - lUMr)Km{Mr)) . 

^ — ' TT ^ — ' 

r7i=— 00 ni=— 00 

Note that the second term in the Green's function (jHUj) can be dropped if the potential 
decreases fast enough at infinity, which is the case here. 

We use the following sum rule for Bessel functions: Ylm=-oo ^m{Fr)Km{Fr') = KQ{F{r — 
r')), therefore, we rewrite the previous expression with two different radii: 



^ 00 

-lim V {Irn{Fr)KUFr') - Im{Mr)K^{Mr')) 



TT r'-+r 

m=— 00 



= - lim[iro(^(^ - r)) - Ko{M{r - r'))]. 

For small r, we have Koi^r) ~ — ln(r) and [KQ{F{r — r')) — Kq^M^t — r'))] = In [y^]- The 
limit r' ^ r is trivial and we get for the whole counterterm: 

h{r)rdr, (A17) 







_F\ 


I 



In 

27r 

which is precisely the counterterm in the Pauli-Villars scheme ()A2|1 

29 



APPENDIX B: ONE-LOOP DIVERGENCES IN PARTIAL WAVES 



The divergent diagrams studied in the framework of Pauh-Villars regularization, see 
Appendix A. 2, A. 3, A. 4, can be recalculated with partial waves for a constant background. 
Their sum is expressed in equation (|37|) and we perform the integration in the following. We 
have: 



(Bl) 



which may be simplified using asymptotic expansions in order for Bessel functions. As 
the divergences are coming from large m, this approximation takes care of the necessary 
contributions: 



Km{Fr) 



1 



2^(^2 + ^2^2)1/4 
2 1 



exp 



exp 



\/m? + F^r^ — m arcsinh ( -z;- ) 



-ym? + Fh^ + m arcsinh 



\FrJ\ 

m 



FrJ\ 



(B2) 



The second term in the propagator ()B1|) is very small if i? ^ 1 and can be neglected. If the 
background is supposed to be constant, it can be taken out of the integral, (j?fj) becomes 



m=—L 
L 

[f{\<P?-v^) + ^-e,,^d,A^) J2 — (Vm2 + FV2 - , 



dr 



m=—L 



(B3) 



where the sum can be converted to an integral: 



2°r / 4L2 \ 



/i/ -1 
— (^Vm^ + F^r^-m]dm 

rewriting (/^ (|0|2 — y"^) + ^e^^yd^A^) as an integral over space lead to 



(B4) 



1. Photon mass term in partial wave 

Finally we recalculate the photon propagator in partial waves. The fermionic contribution 
to the photon propagator comes from three diagrams. The first one reads 



j-R 

V / 27rrrfrG"(r,r)tr(]l)Aj(r). (B5) 
™ Jo 
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This integration is precisely the same as (|H7|). and the result is: 

^2 rR 



e 



l + log(^) 



(B6) 



The second diagram is 




with vertices ie'j^A^d^. In our case, A^(9^ 



Ardr + Ag^dg with A^ = and ^dg is replaced by y for the partial wave m. We further 
assume that Ag is constant over all space. The above diagram gives 



leM^DY, J dhdh'Ag'^G^{r,r')Ag'^G^{r',r), 



(B7) 



where G"^{r,r') is given by ()B1|) . We are interested in large m contributions, and therefore 
we use the asymptotic formulas ()B2l IB2|) for Bessel functions in the propagator. We are 
also interested in the limit R ^ 1, therefore we drop once again the second term in the 
propagator. After some calculations we get : 

11 1 

^^^'^^ = A(9^^2 I 2 . ^1 1 r^T^ ^Tl^ ^^^P [^9{m, r, r')] 9{r' - r) 

Qy.^\lg{m,r\r)\Q{r - r')} . (B8) 

with 

V VP? + F^r^ — V + F^r'^ — m arcsinh V m arcsinh ) . 

Fr Fr' / 

The dominant contribution comes from diagrams with r = r'. Expanding in powers of r — r' 
and performing the integrations, we get for ()B7|) : 



m 

+ / (ir'exp 
Jo 

2 L 



^ m^dr 
+ F^r^ 



R 



dr' exp 



Vm2 + F2 



IY''2i ( ly^ 



m=—L 



1 



— A« 2 / m' 



-A 



FVm2 + ^2/^2 
1 1 



F2m F^y/WTFm^ 



dm 



4 



-1 + lot 



e r .2^2 1 

4 .L ^ 27r 



\F^R^ ) 



+ 0{L- 



(B9) 
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The third diagram is: V/\/V with vertices —ief(j)'jf^Afi. 




y2 



tr(72)Aj02 / rdrr'dr'{27TY{G"'{r,r')Y. (BIO) 



Using an asymptotic expression for the propagator ()B8|) as before, and doing the integration 
in a similar way, we get: 



e 



2 /--R 



/,2 



^ „ , , A^rfV. (Bll) 

The very same way we can recalculate the diagrams ()A12l IA14|) to find respectively 
4^ /cf \9ti4>\'^d'^r and Af^{(f)*df^(f) — (j)d^(f)*)d'^r. These three last expressions can be 

rewritten into a covariant derivative -^^\{df^ — ieA^)(j)\'^. The first two do not cancel com- 
pletely and a term ^A'j^ needs to be subtracted from the action, to get a gauge invariant 
regularization (0^. After this the physical vector boson mass is given by ()A16j) . 



APPENDIX C: EXCHANGING THE LIMITS 



Two limits were considered in the determinant calculation, the limit of infinite volume 
{R oo) and the limit of infinite cutoff in the sum over the partial waves {L oo). The 
order of limits specified in equation ()42p . that is to say take L oo first and then R — > cxd, is 
essential. In this Appendix, we calculate the determinant in the case of vanishing instanton 
core size^° and consider what would happen we commute the limits. In this simple case 
everything can be done analytically; the counterterm to the scalar field mass vanishes, 
because of zero core size. The result for the sum over non-zero partial wave m should be 
finite after removing counterterms related to vector fields. The solutions in the case n = —1 
with boundary conditions ()44|) are: 

(r) = Im-i/2{Fr) ^^^^^^ , {r) = Im{Fr), 

(r) = W/2(Fr) j.^^^^^ , (r) = /„(Fr). (CI) 

Taking a zero instanton core size lead to nornialization problem for the zero-mode. This is not essential 
for our purposes, and it is possible to reproduce all these calculations more rigorously considering a "step" 
core, where /(r) = A{r) — 0, r < S; and then consider S ^ 0. However the calculations are tedious and 
the same conclusions remain. 
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Using 

the determinant for i? — oo is given by: 

...^m. _ vl^r""(oo)^r""(oo) _ r(m + l/2)r(m + 3/2) 

W^' (oojw^' (oo) T[m + lY ^ ' 

= 1 + h C(m ^). 

4m 

Clearly (l + ^) diverges. Note that using this method for the complete numerical 
calculation, the very same divergence remains after removing the ultraviolet counterterms. 

With the second method, using asymptotic expansion ()B2j) for large m and finite radius, 
we get: 

det[M-] nm+l/mm + 3/2) ^^^-.^ + ^^^.^-^ 

^ ^ r(m + l)2 \ Am ^ ' ^ 

= l + 0(m-2). 

which gives a convergent product. This shows that, also in this simple case, we have to 
perform the sum over m to infinity before taking R ^ oo, otherwise we do not get a sensible 
answer. 



APPENDIX D: DETERMINANT AT SMALL FERMION MASS 

Observation of numerical results shows a power law behavior of the determinant for small 
fermion mass. More precisely this power law comes from the partial determinant a/ det'M^, 
where we remove the zero-mode. It is also this contribution that provides the dimension 
of mass~^ for the determinant. It would be interesting to find this behavior by analytical 
calculations. To this end we will use another method H than ^ to remove the zero 
eigenvalue. 

The zero mode wave function which vanishes at the boundary is noted "^oir) and ^o{r) 
shall be the other solution of the second order differential equation ()55|1 : 



^o(r) = e--^o'^'^'f("'), with g{r)=Ff{r) + -A{r), (Dl) 
. L r' 
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This last solution is not normalizable and the constant a which defines the integral is arbi- 
trary. We consider the system to be in a spherical box of radius R. The actual solution, 
which vanishes at the boundary, is not '^o{r) anymore but ^E'a(^); which has a non-zero 
eigenvalue. '^\{r) can be found with the help of perturbation theory: 

^A(r) = ^o(r) - A / tcit[$o(r)^oW-^o(r)$o(t)]^o(t), (D3) 

where the two solutions ()D1UD2|) are normalized so that their Wronskian is 1/r exactly. We 
replace '^\{R) = in the previous equation, this yields 

^ = ^rTWv ^(^)= r tdt[^,{R)^o{t) - ^o{R)^o{t)]'^^{t). 

Jo 

Then the determinant with lowest eigenvalue omitted is 

^o(^) 1 h{R) 



det\M° 



In order to find an analytical approximation for this last expression, we use the following 
approximate profile for the instanton: 

e^r, r < 1/e, 
1/r, r > 1/e, 



eA{r) 



eFr, r < 1/e, 

Ff{r) = { ' - / ' (D4) 
F, r > 1/e. 

Note that the powers of e are introduced for dimensional reasons, the asymptotic behavior 
is exact and the behavior near the center is closely resembling the instanton core. The 
solutions (IDll ID2I1 become 



exp (-ie(e + 2F)r2) , r < 1/e, 
^exp(-i + £-Fr), r>l/e, 

iexp (-ie(e + 2F)r2) 
$o(r) = <j X [Ei (-|e(e + 2Fy) - Ei (-^)] , r < 1/e, 
^ exp(-Fr+i-ii) ^g^p^2Fr) - exp(2F/e)) , r > 1/e. 

Using asymptotic expansions and neglecting parts decreasing as exp(— Fr), the primed de- 
terminant yields: 

det'(M?) = yfexp(i-^)j^\SW** 
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That is to say, for dimensionless variables: 

F^detMl ~ 0.805 (^^^ . (D6) 

It can be compared to numerical results for the partial determinant detM° with which it 
agrees to few percents. The discrepancy comes from the approximate estimate ()D4|) done 
for the instanton profile. The power law behavior is confirmed in figure Note that the 
constant y^ '^^pC^/^) ^ 0.805 in ()D6|) is not expected to match the constant found in the fit 
of figure 5, where the complete determinant was plotted. 
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